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ABSTRACT 



£C) , Aims. We investigate the impact on Fe abundance determination of including magnetic flux in series of 3D radiation-MHD simula- 

tions of solar convection which we used to synthesize spectral intensity profiles corresponding to disc centre. 

Methods. A differential approach is used to quantify the changes in theoretical equivalent width of a set of 28 iron spectral lines 
spanning a wide range in wavelength, excitation potential, oscillator strength, Lande factor, and formation height. The lines were 
P*l \ computed in LTE using the spectral synthesis code LILIA. We used input magnetoconvection snapshots covering 50 minutes of solar 

■ evolution and belonging to series having an average vertical magnetic flux density of (B vert ) = 0, 50, 100, and 200 G. For the relevant 

calculations we used the Copenhagen Stagger code. 

Results. The presence of magnetic fields causes both a direct (Zeeman-broadening) effect on spectral lines with non-zero Lande factor 
and an indirect effect on temperature-sensitive lines via a change in the photospheric T-r stratification. The corresponding correction 
^ ■ in the estimated atomic abundance ranges from a few hundredths of a dex up to |Alog e(Fe) G | ~ 0.15 dex, depending on the spectral 

j.^ ■ line and on the amount of average magnetic flux within the range of values we considered. The Zeeman-broadening effect gains 

relatively more importance in the IR. The largest modification to previous solar abundance determinations based on visible spectral 
lines is instead due to the indirect effect, i.e., the line-weakening caused by a warmer stratification as seen on an optical depth scale. 
Our results indicate that the average solar iron abundance obtained when using magnetoconvection models can be ~ 0.03 - 0. 1 1 dex 
higher than when using the simpler HD convection approach. 

Conclusions. We demonstrate that accounting for magnetic flux is important in state-of-the-art solar photospheric abundance deter- 
minations based on 3D convection simulations. 



Key words, magnetohydrodynamics (MHD) - radiative transfer - line: formation - Sun: abundances - Sun: granulation - Sun: 
photosphere 



1. Introduction line formation is of secondary importance. The space- and time- 

' averaged solar spectrum commonly used for abundance determi- 

^ , The chemical composition of stars is derived through the anal- nationfj wm ^ fect fee wei hted towards granules , owin g to the 

. ysis of their emitted electromagnetic radiation. Thanks to ap- , surface ^ ^ ^ coyer ^ a tQ int ular 

m > ■ propnate modelling, the presence of lines of a given chemical lm where netic field concen trations tend to be found. The 

^ . element in the observed spectra can be interpreted in terms of idea . g ^ for Ms Kason ^ { Qf etic flux over the 

1 a b unda " ce "> the object. The line strength as quan- a e ent ctmm (an(f in particulal - abundances) is 

d ■ tified, say, by the equivalent width, needs to be matched by the ^ minimized . According to this assumption, the uncertainty 

corresponding theoretical prediction A more detailed approach introduced b usi magnet ic-free models would be less impor- 

invo ves matching the line profiles themselves instead ^ than Qther more obyious SQurces of such ag atomic 

As our closest star, the Sun allows the most detailed obser- ^ uncertainties related t0 obse rvational data, treatment of 

vations of the processes taking place in an object of its type. non _ LTE effectSi and so on . MoreoV er, the assumption that mag- 

The atmospheric models used to describe the physical condi- netic fieWs couW fee disre ded a i toget her in abundance studies 

tions of the solar plasma have become more and more real- hag been co nded with the need of s i m piif ying the simula- 

istic through the years, in particular, accounting for 3D mho- don ^ thesis computations and subsequent anal- 

mogeneities and for the need to have good granulation statis- ig Because of tQ Ms situation; when advancing from the i D 

tics. Still, magnetic fields are usually neglected in solar abun- lane - allel approx i mat i on to a 3D approach, non-magnetic 

dance studies, on the ass umption that their effect on spectral model atmospheres were used . Nevertheless, two major argu- 

* ai ffiv , j ^ , r c, i t-,i . c .. ments can be put forth against that simplification. 
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(STARPLAN), K0benhavns Universitet, 0ster Voldgade 5-7, DK-1350 1. From an observational point of view, evidence is accu- 

K0benhavn 0, Denmark mulating about a large amount of hidden magnetic flux being 
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present in the Sun, even in "quiet" regions (see Sanchez Almeida 
& Martinez Gonzalez 201 llfor a recent review). It is now known 
that weak magnetic fields cover most of the Sun's surface, with 
indications th at a flux density of order 10 2 G may be a good ref - 
erence value dTruiillo Bueno et all 20041; INordlund et all 20091). 
In regions with strong solar magnetic fields (e.g., in plage re- 
gions), the fraction of magnetic field concentrations with of the 
order of 1 kG in intergranular lan es tends to increase (see e.g. 
Fig. 46 in INordlund et al. 2009). In observed plage regions, 
the granulation presents a lot of fine structuring at high spa- 
tial resolution: isolated bright points, strings of bright points and 
dark micro-pores, "ribbons", or more circular "flower" structure s 
(iTitle et al.l 1 991 iBerger et al.l2004tlNaravan & Scharmerl2010h . 
Granules become smaller than in th e quiet areas jMuUeiHl977t 
ISchmidtetalIl988l : lTitle etalJ 19921) . and intergranular lanes are 
characterized by the presence of micro-pores. In agreement with 
observations, numerical 3D convection models predict that gran- 
ulation patter n is profoundly affected in regions with high flux 
densi ty (e. g. ICattaneo et all 2003I IVogleri 2005, Ch eung et all 
2008 ), with reduced size of the upflowing granules and widened 
intergranular lanes. In the extreme case of sunspot umbrae, the 
convection is s eriously impaired and take s the form of narrow 
upflow plumes dSchiissler & Vogle r 2006). 

2. The reasoning given for why HD should be a good ap- 
proximation for line formation in the Sun neglects the fact that 
the photospheric material becomes more transparent in mag- 
netic concentrations owing to their lower density, allowing one 
to probe deeper, i.e., into hotter layers of the atmosphere (Stein 
& Nordlund 20001) where the flux tubes are also heated through 
heat influx from the surrounding material dSpruitl 19761). Weak 
spectral lines will then experience brightening of their core. The 
we akening of neutral iron spectr al lines was, for example, used 
by ISchuessler & Solankil ( 1988) to predict a higher continuum 
intensity for magnetic elements. Observed photospheric bright 
points may then be identified with magnetic flux concentrations. 
Magnetic fields can therefore act on temperature-sensitive spec- 
tral lines not only directly, i.e., via Zeeman broadening, but also 
via an indirect effect associated with the temperature stratifica- 
tion in the line-forming regions. Only a comprehensive study us- 
ing realistic 3D magnetohydrodynamical simulations can there- 
fore confirm whether the effect of magnetic fields on spectral 
line formation and abundance derivation is truly negligible. 

Including magnetic fields is clearly an important new step 
in the long-term ev olution of abund ance studies. Following pi- 
oneering work (e.g. INordlund 1984), 3D radiative magnetocon- 
vection simulations that explore realistic configurations and pro- 
vide an excellent match to the photometric and spectroscopic ob- 
se rvations were produced during the past decade, as summarized 
by INordlund et al.l (2009) The resulting 3D model atmospheres 
(albeit with zero magnetic field) have been used in abundance 
determination work to realistically account for the departures 
from homogeneity in line formation. Along with the use of better 
input atomic data, these advances led to remarkable agreement 
between the observational data and the theoret ically predicted 
shapes and asymmetries of the spectral lines (lAsplund et alJ 
l2000al : IShchukina & Truiillo Buenol 2001l). but only if decreased 
abundance of several chemical elements in the Sun was assumed 
(lAsplund et alj 2000al; lAsplundl 2005l; lAsplund et al.l 20091) com- 
pared to the previously commonly adopted standard reference 
values inferred via a ID static, homogeneous, and plane-parallel 
LTE analysis. The 3D-derived abundances give smaller line-to- 
line scatter compared to the ID results at the same time as ren- 
dering the microturbulence and macroturbulence free parameters 



unnece ssary. Better agree ment with solar neighbourhood abun- 
dances dLodders et al. 20091) is also obtained. 

While all these results are very encouraging, pressing out- 
standing issues remain. Of particular importance among them 
is the neglect of magnetic fields in the convection models men- 
tioned above, an issue that we approach in this paper. The role of 
magnetic fields in abundance determin ations has been initially 
investigated in ID model atmospheres dBorrerol 2008) finding 
non-negligible effects, but it was not until very recently (Fabbian 
et al. 2010) that an investigation set out to explore the effects 
in the 3D MHD case. The results of that exploratory paper in- 
dicated that, even in the magnetically quiet Sun, the neglect 
of magnetic fields may question the few-centidexes accuracy 
claimed in previous works. 

In the present paper, we achieve the first comprehensive 3D 
study of the solar iron abundance that includes magnetic fields. 
We focus on iron because of the large number of spectral lines 
for this element in the visible and IR with different direct mag- 
netic sensitivity (i.e., differing Lande factors gi) and also be- 
cause the solar Fe abundance is the basis for most abundance 
work in the astrophysics literature. Iron is in fact of central his- 
torical importance as a proxy for the total metal content of a 
given stellar object and as the reference element in the elemental 
abundance ratios used in studies of galactic chemical evolution. 
The Fe solar abundance was clai med in the 80's to lie anywhere 
between loge(F e) M = 7.70 d ex (Blac kwell et~ai] 19841, 19861) 
and ~ 7.48 dex dHolweger et al.l 19901) The use of 3D model at- 
mospheres and improved observational, atomic, and modelling 
data led to a drastic reduction of the large line-to-line abun- 
danc e scatter, whereby a low abundance of log c(Fe) n = 7.45 
dex (lAsplund et a l. 2000b) was det e rmine d, consistent with 
the meteoritic value of lLodders et al.1 (2009) The question re- 
mains whether that mean value can withstand the inclusion of 
magnetic effects in the 3D model atmospheres by using mag- 
netoconvection (as opposed to purely HD) models. 

In this paper we find confirmation of the importance of mag- 
netic fields in determining the solar Fe abundance. For the F e 
spectral lines we have in common with lAsplund et all (2000b) 
we find abundance corrections of Alog e(Fe) = 0.05 - 0.09 dex 
when considering magnetoconvection models with an average 
unsigned magnetic flux of (B vert ) = 100 G. 

The neglect of magnetic fields is not the only remaining 
problem in the abundance determination research field. Other 
open questi ons are (i) the friction with helioseismic inversion 
constraints (lAntia & Basu 2005 ; Delaha ve & Pinsonne ault 2006 ) 
caused by the up to ~ 30% smaller estimate currently widely 
adopted for t he solar photospheri c metal content; (ii) several 
authors (e .g., Caff au et al] 20071; lAv res 2007 ; ICaffau et al.l 
20081; lAvresI 20081; ICaffau et al.l 20101; ICaffau et al.l 201 II) have 
suggested that the revision in the abundance of a number of 
chemical elements may be too drastic, for example due to 
slightly too steep temperature gr adients in the employed mod- 
els' continuum-forming layers dTruiillo Bueno & Shchukinal 
2009 ) or too-low tem peratures in th eir middle photospheric lay- 
ers (lAvresI 20121); (iii) ICaffau et al.l (20091) claim that the revised 
solar abundances are likely mostly the result of differences in 
equivalent width and/or line profile fitting, as well as in atomic 
data and blends, rather than granulation effects and hydrodynam- 
ical modelling. 

The layout of the paper is as follows. Section [2] deals with 
the details of the atmospheric models employed in this study 
and with the method used in computing theoretical spectral line 
profiles. Section|3]presents the abundance correction results that 
we find for the different Fe lines studied here. In Sect. |4]we draw 
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consequences specifically with respect to the solar iron content. 
Section [5] discusses the more general implications in terms of 
derivation of the elemental composition of stars, and Sect.|6]lists 
our conclusions. 



2. Atmospheric models and computation of 
spectral lines 

2.1. Method 

We used simulations covering several hours of solar time that we 
obtained by usi ng the 3D radiative-hydrodyna mics Copenhagen 
STAGG ER code (|Gal sgaard & Nordlund 1996 also, see discus- 
sion in Fabbian et al. 2010 and the recent description in Beeck 



et al. 20121). The simulation box covers 6.0x6. OMm 2 horizon- 
tally and extends in height for -2.8 Mm. These datacubes are di- 
vided into 252x126x252 cells, with non-uniform spacing in the 
vertical direction, yielding a resolution in height of ^ 15 km at 
the photosphere. In the vertical direction, after avoiding the ver- 
tical ghost zones necessary in the computations, only ~ 2.5 Mm 
actually contain valuable physics information. Horizontally pe- 
riodic boundary conditions are imposed for all variables on the 
sides of the simulation box. 

The initial setup configuration for the magneto-convection 
simulations is a uniform vertical magnetic field. While diverg- 
ing convective upflows will tend to sweep and c oncentrate mag- 
netic fields to intergranular downflow lanes (e.g.. lNordlund et all 
2009), the horizontal average of the vertical component of B 
must remain equal at all depths and times throughout the sim- 
ulation to the initial value chosen for it in the relevant series, 
< B vert >= B f= o = const, and thus characterizes the different 
magnetic series in this study. We subdivide our results into se- 
ries for initial vertical field of G (the HD series) and 50, 100, 
and 200 G for the MHD series. For each series we wait until 
a statistically stationary state has been reached, which typically 
takes some 0.5 solar hours after introduction of the initial ver- 
tical magnetic field. We then store snapshots of the subsequent 
temporal evolution at intervals of 30 s solar time. 

F or the a-posteriori spe ctral synthesis we used the LILIA 
code dSocas-Navarr o 2001 ). A subdomain of the computed box 
extending from 425 km above to 475 km below the average 
i"5oonm = 1 level was used for the synthesis; the spatial resolu- 
tion in that domain was increased by interpolating to a uniformly 
spaced grid with resolution ~ 7.8 km. In Fig. Q] we show the 
MHD-HD differences in the spatially and temporally averaged 
temperature stratification. The average temperature stratification 
is different for the different series, with a tendency in layers with 
logT5oonm < to an increasingly large departure from the HD 
temperature stratification as the initial magnetic field becomes 
higher, following the expected behaviour (see Sect.[TJ. We stress 
that this effect appears due to the changed opacity caused by the 
decreased pressure (lower density) inside strong magnetic flux 
concentrations (magnetic plus gas pressure being nearly equal 
to the gas pressure in the surrounding unmagnetized plasma, e.g. 
iNordlund et al. 2009 ). The photospheric temperature increase in 
MHD is therefore not related to the magnetic energy build-up by 
photospheric motions, which has been hypothesized to heat the 
corona after local release via magnetic reconnection. 

The synthetic spectra were obtained only for solar disc cen- 
tre. After verifying that the influence on the correspondingly de- 
rived theoretical equivalent widths was negligible, for the spec- 
tral synthesis we decreased the horizontal resolution by taking 
one out of every four grid points. We selected a subset of snap- 
shots covering the last 50 minutes of the statistically stationary 
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Fig. 1: MHD-HD difference of the temperature stratification obtained 
for layers close to photospheric line-forming regions after averaging 
horizontally on surfaces of equal optical depth and temporally over the 
(63x63 column) input snapshots used in the spectral synthesis. Blue, 
green and red lines represent the difference in temperature with respect 
to the HD case, against log 10 r^oonm) (bottom horizontal axis), for the 50 
G, 100 G and 200 G cases, respectively. The top horizontal axis shows 
the corresponding physical depth in Mm in the HD case. 



regime of convection, separated by a 150 solar seconds time in- 
terval. We therefore keep 21 snapshots per series. 

Spatially- and temporally-averaged synthetic spectra were fi- 
nally computed and then employed for a differential analysis of 
the impact of different amounts of magnetic flux on abundance 
determinations. 



2.2. Line parameters 

The 28 iron lines we chose for the spectral synthesis are listed 
(with their parameters) in Table Q] The spectral lines cover a 
wide range in rest wavelength {A), with 23 of them in the vis- 
ible and 5 of them in the near-IR spectral region. To detect 
possible dependences of the final abundance corrections on the 
atomic parameters, the spectral lines were chosen so as to in- 
clude a wide range of values of lower transition level excitation 
potential (0.1 < xi < 6.5 eV), of transition oscillator strength 
(-5.0 < log gf < 0.4), and of Lande factor (0 < g L < 3). The 
values adopted are from the VALD (Version 0.4) and NIST ASD 
(Version 4) database^ 

In Table Q] we assign different identifiers for subgroups of 
lines, to reflect the organisation that will be used in the dis- 
cussion of results (Sect. |3j. We make a division into Zeeman- 
insensitive (gL = 0) lines (Sect. 13.2b - identified in the table 
with the label "0L" - and in lines with non-zero Lande factor 
(Sect. [33}. The latter are additionally subdivided into lines with 
high values of the Lande factor ("LL" in Table [p, i. e., having 
gi. > 1 (see Sect. 13.3.2b . lines in common witl uAsplund et all 
(2000 b|) (Se ct. [3331 ) - marked in the table with "A" - and IR 
lines (13.3.4b . identified with "IR". The spectral features listed in 



The Vienna Atomic Line Database (VALD) is available as an on- 
line interface at http://vald.astro.univie.ac.at/~vald/php/vald.php; 
the National Institute of Standards and Technology 
Atomic Spectra Database (NIST ASD) is accessible at 
http://www.nist.gov/pml/data/asd.cfm 
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Table 1 : Atomic parameters for the 28 iron spectral lines included in the 
current study. 



A 


Xi 


log,?/ 


gL 


ID 


[nm] 


[eV] 








410.091 


2.4530 


-4.244 


0.0 


0L 


415.025 


3.4300 


-1.260 


0.0 


0L 


423.027 


4.2310 


-2.440 


0.0 


0L 


430.603 


4.2180 


-1.922 


0.0 


0L 


437.113 


4.2600 


-3.306 


0.0 


0L 


486.364 


3.4300 


-1.663 


0.0 


0L 


524.705 


0.0873 


-4.946 


2.0 


A 


525.021 


0.1213 


-4.938 


3.0 


A 


525.065 


2.1979 


-2.181 


1.5 


LL 


543.452 


1.0110 


-2.122 


0.0 


0L 


557.609 


3.4302 


-0.940 


0.0 


0L 


579.119 


5.0200 


-2.067 


0.0 


0L 


585.923 


4.3010 


-3.076 


0.0 


0L 


594.036 


5.0100 


-2.260 


0.0 


0L 


601.392 


5.1000 


-2.133 


0.0 


0L 


608.271 


2.2230 


-3.572 


2.0 


A 


611.332 


3.2210 


-4.230 


0.57 


11 


616.946 


3.4170 


-4.009 


0.0 


OL 


617.334 


2.2227 


-2.880 


2.5 


A 




9 99^0 


J.J7U 


i n 
1 ,u 


A 
r\ 


630.150 


3.6540 


-0.718 


1.67 


LL 


630.249 


3.6860 


-0.973 


2.5 


LL 


633.084 


4.7330 


-1.740 


1.23 


LL 


1558.826 


6.4740 


0.391 


1.5 


1R 


1564.851 


5.4260 


-0.670 


2.98 


1R 


1565.287 


6.2461 


-0.170 


1.55 


1R 


1566.202 


5.8280 


0.382 


1.47 


1R 


1566.524 


5.9790 


-0.338 


0.67 


IR 



Notes. The columns list, respectively, the rest wavelength, lower 
transition level excitation potential, oscillator strength, Lande factor 
of the iron spectral lines studied here, and the identifiers employed in 
Sect.[3]for a useful subdivision into different groups of lines. 



Table[T]are all lines of neutral iron, except for the singly-ionized 
61 1.332 nm iron line, marked in the table with "II". The group- 
ing just discussed will help to more clearly explain the direct and 
indirect effects of the magnetic flux on spectral lines. 

The iron lines chosen for this study probe quite different 
depth ranges in the photospheric snapshots, thus jointly cover- 
ing a domain from the high layers in the simulation box to as 
deep as around 5 km above log 10 T^soonm) = 0. Their approx- 
imate heights o f formation (e. g. see values in Gurtovenko & 
Kostik 1989, Table 2) make them generally safe for this study 
since they are supposed to be forming inside the simulation do- 
main. 

The collisional Van der Waals broadening/da mping formula 
was used following the classic approximation of lUnsoldl ( 19551) 
with no enhancement. Our tests revealed that, while the exact 
choice of the enhancement factor will influence the resulting 
equivalent width and is thus an important choice for line profile 
fitting and absolute abundance matching, it does not affect a dif- 
ferential abundance investigation, like ours, to a significant level. 
For instance, the abundance correction derived from the compar- 
ison of the HD and MHD equivalent widths becomes less severe 
by less than 0.01 dex in absolute value when the enhancement 
factor is changed from 1.0 to 2.5. 

The lines were synthesized at a spectral resolution of 0.50 
pm. We found this to be a good compromise for achieving high 
spectral resolution at an acceptable computational effort with a 
resulting equivalent width accuracy better than 0.001 pm. For 
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Fig. 2: Top panel: comparison of our computed solar disc centre contin- 
uum intensity in the non-magnetic and 200 G case (solid black and red 
curves, respectively) to data from the literature. The long-dashed curve 
is based on the results obtained bv lTruiillo Bueno & Sh chukina (2009) 

w lAs 



using a single 3D HD snaps hot bvlAs plund et al. (2000a), the dotted 
curve with circles is based on lAvres et a l. ( 2006) , the short-dashed curve 
refers to the semi-empirical ID MACKKL mod el of the "quiet" Sun at- 
mosphere (see Table II of iMaltbv e t al. 1986); fin ally, b lack crosses 
indicate the observational data from lNeckel & La bs ( 1984, Table VII) 
Bottom panel: intensity difference (in percentage) between our results 
and the observational data (black crosses=HD-obs; red crosses=200 G- 
obs). 



each spectral line we used as many spectral points around the 
central wavelength as needed to reach convergence at least at 
the 10~ 2 level when calculating the equivalent width. For the 
strongest lines this implied using up to 1000 wavelength points 
(or ~ 250 pm on each side of the rest wavelength). 

2.3. Continuum intensity 

Figure [2] (top panel) shows the level of agreement between our 
calculated continuum intensity data (based on a space and time 
average for the same snapshots of the HD series that were em- 
ployed in the spectral synthesis) and different theoretical and 
measured values, all being spatially-averaged means at solar 
disc centre (ji — 1) and given in absolute units. Our continuum 
intensity results for the HD and 200 G series are seen to lie be- 
tween the 3D HD values that lTrujillo Bueno & Shchuki na ( 2009) 
derived based on a single snapshot from lAsplund et al.1 (2000a) 
and the va lues from the solar photospheric "thermal profiling" 
analysis of lAvres et al.l (2006L see their Table 1) Our model lies 
in a band of < 5% below the semiemrjirical model of Maltby 
et al. ( 19861) known as MACKKL (short-dashed curve in Fig.|2]i, 
which was purposely built at the time precisely from continuum 
observations. Therefore, that our results are, in particular, close 
to the latter, confirms the goodness of our 3D MHD simulations. 
They also agree well with the observational data from Neckel 
& Labs (1 984. Table VII) An d with the continuum distribution 
of lShchukina & Sukhorukovl (2012l) (not shown in our figure to 
avoid the plot becoming overcrowded) , which is based on the lo- 
cal dynamo simulation of lVogler et al.l ( 20051) The bottom panel 
of Fig.|2]then shows the intensity difference (in percentage) be- 
tween our results and the observational dataset of lNeckel & Labsl 
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( 1984) The cases shown are our HD and 200 G series (the differ- 
ence with observations being marked in the plot as black and red 
crosses, respectively). The time-averaged HD and MHD mod- 
els both show a good fit to observations. The uncertainties in- 
volved in a correctly computed continuum are of the same order 
as the differences we find with respect to available literature con- 
straints. 

3. Abundance correction results 

3.1. Overview 

The procedure used to obtain abundance corrections is as fol- 
lows. From the (M)HD snapshots we obtained synthetic spec- 
tra for each of the columns in the snapshots and for all of the 
lines in Table [TJ using a standard reference abundance value of 
log e(Fe) = 7.50 dex (although the precise value is not crucial 
here, given the differential nature of our study). Then, we cal- 
culated space and time averages (i.e., horizontally in each snap- 
shot and for all snapshots) separately for each of the series (HD, 
50 G, 100 G, 200 G). We then computed the equivalent width 
of the average profiles and compared them across the series. 
Additionally, to determine abundance corrections, we ran spec- 
tral synthesis calculations for the HD case with an input Fe abun- 
dance changed by -0.10, -0.06, -0.02, +0.02, and + 0.10 dex 
compared to the reference value mentioned above. The precise 
abundance correction to use in the HD series to match the equiv- 
alent width of the profiles for the MHD series was then found by 
interpolation between the values obtained for these fixed abun- 
dance jumps. In fact, the abundance correction to adopt as a re- 
sult is actually the negative of the value obtained through the 
interpolation. The reason for the change of sign is that if, e.g., a 
given spectral line tends to weaken for non-zero magnetic flux, 
then when neglecting magnetic fields the equivalent width mea- 
sured from observations will be matched with a lower (i. e. un- 
derestimated) elemental abundance. 

Finally, to separate the magnetic broadening (Zeeman split- 
ting) effects from the influence of the changed temperature struc- 
ture in the MHD models, we ran further LILIA spectral synthesis 
calculations based on the MHD models (in particular, keeping 
the temperature and density stratification in them) but artificially 
setting B — 0. We expect that the line formation in this case 
will contain the indirect effect of the magnetically-induced tem- 
perature stratification change, but of course, without any direct 
Zeeman broadening. 

In the following sections, we present the results of this study 
by dividing the spectral features in our line list into two broad 
groups. The first group (Sect. l3~2l contains the 13 lines with no 
Zeeman sensitivity (gL ~ 0), while the second group (Sect. 13.3b 
is composed of the 15 lines with at least some direct sensitivity 
(0.5 < gh < 3.0) to the magnetic field. 



3.2. Zeeman-insensitive spectral lines (g L = 0) 

Spectral lines with gi — are indicated with "0L" in Tables [TJ 
and |2 Any magnetic effects on these lines must necessarily 
come from the change of temperature and density stratification 
alone. The magnitude of the effect then depends on the forma- 
tion height and temperature sensitivity of the line. The behaviour 
of the calculated equivalent width with respect to magnetic flux 
is illustrated for two lines of this group in Fig.[3](left and middle 
panels). The corresponding abundance corrections are included 
in Table [2] together with the predicted HD equivalent width for 
each spectral line (which can provide an idea of their approxi- 



mate for mation region, since the formation height is not avail- 
able from iGurtovenko & Kos tik 1989 for all lines in our list). 
The abundance corrections turn out to be small +0.02 dex) 
for Fei 423.027 nm (left panel), but reach up to ~ +0.10 dex 
for Fe 1 543.452 nm (middle panel) in the 200 G case. The hotter 
temperatures that the radiation "feels" in the MHD models make 
(for a fixed adopted abundance) these spectral lines increasingly 
weaker the higher the unsigned flux of the different MHD series, 
in agreement with their temperature sensitivity. By checking the 
position of the dashed horizontal lines, we see that the abun- 
dance of the HD case must be decreased appreciably for it to 
match the equivalent width found for the MHD cases. Using the 
non-magnetic case to match measured equivalent widths from 
observed spectra, one would then derive a low abundance esti- 
mate. The latter will in fact allow to match predicted and mea- 
sured equivalent widths because the spectral line is intrinsically 
stronger in HD. This finally translates into a positive abundance 
correction to the purely HD case, because the abundance values 
obtained for it must be revised upward in order to match ob- 
servations while also endeavouring to account (precisely via the 
use of the relevant abundance correction) for the missing effects 
(instead naturally present in MHD) due to magnetic fields. 

The other Zeeman-insensitive spectral lines in our sample all 
have resulting abundance corrections within the above values. 
For example, they reach up to only ~ +0.03 dex for the Fei 
585.923 nm feature, but up to ~ +0.09 dex for the 486.364 nm 
line. 

As seen in Fig. [3] (left and middle panels), the Zeeman- 
insensitivity of all absorption lines in this category is confirmed 
by the absence of any equivalent width variation when artificially 
setting B — in the spectral synthesis. 

3.3. Spectral lines with g L ^Q 
3.3.1. Fen 611.3 nm line 

Figure [3] (right panel) shows the equivalent width results for 
the only Fe n spectral line included in this study (marked with 
"II" in Tables [TJ and The corresponding corrections to the 
HD-derived abundance are modest. Since Fen lines are weak 
(our predicted HD equivalent width for this line is 0.92 pm, as 
listed in Table |2]l and generally formed rather deep in the photo- 
sphere, in the region where the temperature modifications due to 
the presence of magnetic fields are small (see Fig. [TJi, the cor- 
responding abundance corrections can also be expected to be 
fairly mild. The formation height of the Fen 611 .332 nm lin e 
lies around 1 10km (see Table 2 in IGurtovenko & Kos tik 1989 ). 
This line, moreover, belongs to the visible wavelength range and 
has a small Lande factor (gL = 0.57), so that, as seen in Fig. [3] 
(right panel), the Zeeman broadening is unable to contribute any 
significant line strengthening. The applicable abundance correc- 
tions thus settle at a maximum of ~ +0.03 dex, see Table [2] 
Singly-ionized iron lines are generally less sensitive to tempera- 
ture, to no n-LTE effects and to details of the model atmosphere 
(see, e. g.. lAsplund et al. 20091) Improved atomic data for them 



are now available (e.g., Melendez & Barbuv[ 2009). It will there- 
fore be of future interest to ascertain, based on a larger sample 
of singly-ionized iron lines, whether they may also be generally 
insensitive to the direct and indirect effects of magnetic fields, 
as the results for at least Fen 61 1.332 nm seem to suggest. This 
would likely confirm Fe n as a more reliable indicator than Fe i 
for the purpose of abundance determinations under simplified 
(e.g., LTE, or non-magnetic) approximations, i.e., when a more 
sophisticated analysis is not feasible. 
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Fig. 3: Left and middle panels: equivalent width results for two of the Zeeman-insensitive (g L = 0) spectral lines in our line list. Right panel: 
equivalent width results for the Fen 61 1.332 nm spectral line. In both panels, filled circles mark the equivalent width for the (M)HD cases, with 
average flux density (< B rerl >= B l=0 = const, i. e., equal throughout the evolution to its initial value) given in abscissas. Empty circles (here, 
essentially overlapping the filled ones) represent the equivalent width results obtained when artificially setting B = for the spectral synthesis. 
The solid horizontal line marks the value of the equivalent width for the HD case. The dashed horizontal lines mark the equivalent width value 
obtained when decreasing the iron abundance of the HD model by 0.02, 0.06, and 0. 10 dex. 
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Fig. 4: Equivalent width results for three of the g L > 1 
spectral features in our line list. The symbols are the same as in Fig.[3](but now filled and empty circles are no longer overlapping). 



3.3.2. Lines with large Lande factor g L 



Spectral lines with g£ > 1 are indicated with "LL" in Tables Q] 
and [2] The abundance corrections for these lines are given 
in Table [2] together with their predicted HD equivalent width. 
Figure |4] shows the equivalent width results for three spectral 
features with large g^, namely the 630. 150 nm, 630.249nm, and 
525.064 nm lines. The abundance corrections for these lines are 
of intermediate magnitude, reaching the value ~ +0.05 dex, 
~ +0.06 dex, and ~ +0.08 dex, respectively. Given their non- 
zero Zeeman sensitivity, the equivalent width of all three lines 
decreases when artificially setting B — for the spectral syn- 
thesis, as seen in the figure. However, this effect remains small 
throughout the visible range, and is not sufficient to counter the 
indirect line-weakening temperature effect. 



Some of the lines marked as A or IR in Table Q] also have 
a large Lande factor gL- However, they are discussed as sepa- 
rate subgroups below, because they are helpful for highlighting 
some of the main results of this study, i.e. the possible impact 
on the accepted value of the solar iron abundance and the more 
noticeable Zeeman broadening effects in the IR, respectively. 



3.3.3. Lines in common with the solar iron abundance study 
of Asplund et al. (2000) 

Figure [5] shows the equivalent width results for the spectral fea- 
tures marked with "A" in Tables Q] and [2] We include them in an 
own separate subgroup, given that they are the ones in common 
with a previous well-known solar Fe abundance study (Asplund 
et al. 2000b) Corresponding abundance corrections and HD 
equivalent widths are listed in Table [2] These spectral lines have 
some of the largest Lande factors among the visible lines in our 
sample. However, since all five are in the visible range, Zeeman 
broadening actually has a little influence on them and can do 
only little to alleviate the line weakening tendency, as is apparent 
by comparing the empty and filled circles in the panels. On the 
other hand, their sensitivity to the magnetically induced changes 
in the average temperature stratifications is high, placing these 
spectral features among those in our line list that have the largest 
abundance corrections (see Table|2]i. Namely, for the 200 G case, 
the corresponding correction to the HD-derived abundance is 
between 0.08 dex for the 608.271 nm line to 0.14 dex for the 
524.705 nm one, which is the maximum in the whole sample we 
studied. 

Fei 525.021 n m (one of the two spectral lines in the "good 
magnetic pair" of lStenflol 19731) has the second most significant 
corrections, Aloge(Fe) = 0.12 dex. Apart from having been 
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Fig. 5: Equivalent width results for the spectral lines in common with Asplund et al. (2000b) Same symbols as in Fig. [3] 



used in the solar iron abundance derivation of Asplund et al. 
(2000b) this line has also recently been t argeted by the IMaX 
imaging vector magneto graph instrument dMartmez Pillet et al.l 
201 l HBonet e t al. 2010) onboard the ballon-borne Sun rise tele- 
scope mission dBarthol et all 201 1 ; So lanki et al.l 20101). 



3.3.4. IR lines 

Spectral iron lines belonging to the infrared wavelength range 
are identified with "IR" in TableQ] Table[3]lists their abundance 
corrections and HD equivalent widths. All of the IR lines in- 
cluded in our line list are strong, 7.6 < Whd < 19.9 pm. Figure[6] 
shows the relevant equivalent width results. As seen, these lines 
behave like those in the visible concerning their sensitivity to the 
warmer temperature structure in the MHD models. However, the 
IR iron lines we selected have quite large Lande factors. Coupled 
with the A 2 dependence of Zeeman broadening, this makes them 
prone to strong direct magnetic effects. The result is that, com- 
pared to the HD case, these IR iron lines tend to generally be- 
come stronger in the MHD models. The corresponding abun- 
dance corrections range from negligible to -0.08 dex. 

4. Consequences concerning the solar Fe 
abundance 

Different neutral and/or ionized iron lines should yield a 
consistent estimate for the photospheric solar abundance. Once 
using appropriately weigthed averaging of the abundance values 
inferred from the various spectral lines of a given chemical 
element, little scatter should in fact result around the "best 
estimate". In this section, we provide average values for the 
abundance correction to be applied to different groups of 
lines. In doing so, we just conform to the usual practice in the 



literature whereby averages of abundance values from multiple 
spectral lines have been routinely taken as the representative 
derived solar abundance. Since there is no obvious way to attach 
a different statistical significance to the abundance correction 
obtained for each line, we simply give equal weight to them 
when calculating averages. 

Spectral lines in the visible: as can be seen in Table [2] the 
abundance corrections for the lines in the visible range vary from 
negligible to as high as +0.14 dex. The unweighted mean of the 
corrections for these lines is shown in the row near the bottom 
marked with () signs and it varies between 0.02 and 0.06 dex 
when going from the 50-G to the 200-G case. If we focus on 
the visible Fe lines marked with "A " in our table (i.e., those in 
common with lAsplund et aL (2000b)), we find that they yield 
significantly higher abundance values when MHD modelling is 
used. In detail, out of these five spectral lines, Fe 1 608.271 nm is 
the least affected, with abundance corrections of 0.02-0.08 dex, 
depending on the magnetic flux considered. The most affected is 
Fe 1 524.705 nm, with abundance corrections of 0.04 - 0.14 dex. 

We may apply the average abundance correction values 
(given in the row marked with Q* in Table [2]i to the abundance 
determinations by lAsplund et al.l (2000bl) for those five lines. 
Their average value was loge(Fe) Q = 7.42 ± 0.03. We would 
then obtain an average solar iron abundance estimate in the 
range 7.45 < log e(Fe) < 7.53 dex, depending on the magnetic 
flux assumed. On the other hand, even if we apply abundance 
corrections separately for the different lines following the results 
in Table [2] and calculate the mean value and rms, the value of 
the latter turns out to be close to the one calculated using purely 
HD models. In detail: the rms is 0.03 dex for the 50 G case and 
0.04 dex for the 10 and 200 G cases, henc e very similar to the 
values obtained bv lAsplund et al.l (2000b) Therefore, the good 
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Fig. 6. Equivalent width results for the IR spec- 
tral lines included in our sample. The symbols 
are the same as in Fig. [3] except that now, the 
dashed horizontal lines indicate the equivalent 
width value predicted when increasing the iron 
abundance in the HD model by 0.02 and 0.10 
dex. This translates to negative abundance cor- 
rections of the same magnitude applicable to 
the abundance predicted using the HD model 
(see discussion in text). 



Table 2: Abundance corrections for the Fe spectral lines in the visible range, for the different magnetic flux cases. 



A 


Alog6(Fe) 


Alog e(Fe) 


Alog e(Fe) Q 


W HD 


ID 


[nm] 


[dex] 
(50G-HD) 


[dex] 
(100G-HD) 


[dex] 
(200 G - HD) 


[pm] 




410.091 


+0.01 


+0.03 


+0.04 


0.55 


0L 


415.025 


+0.02 


+0.04 


+0.06 


7.94 


0L 


423.027 


+0.01 


+0.01 


+0.02 


0.69 


0L 


430.603 


+0.01 


+0.02 


+0.03 


1.98 


0L 


437.113 


+0.01 


+0.01 


+0.02 


0.10 


0L 


486.364 


+0.02 


+0.05 


+0.09 


6.60 


0L 


524.705 


+0.04 


+0.09 


+0.14 


6.74 


A 


525.021 


+0.04 


+0.08 


+0.12 


6.72 


A 


525.065 


+0.03 


+0.06 


+0.08 


9.56 


LL 


543.452 


+0.04 


+0.07 


+0.10 


24.72 


0L 


557.609 


+0.02 


+0.05 


+0.07 


10.91 


0L 


579.119 


+0.01 


+0.02 


+0.03 


0.37 


0L 


585.923 


+0.01 


+0.02 


+0.03 


0.18 


0L 


594.036 


+0.01 


+0.02 


+0.03 


0.25 


0L 


601.392 


+0.01 


+0.02 


+0.03 


0.28 


0L 


608.271 


+0.02 


+0.05 


+0.08 


3.84 


A 


611.332 


0.00 


+0.01 


+0.03 


0.92 


11 


616.946 


+0.01 


+0.03 


+0.05 


0.16 


0L 


617.334 


+0.03 


+0.07 


+0.11 


6.98 


A 


624.065 


+0.03 


+0.06 


+0.10 


4.76 


A 


630.150 


+0.02 


+0.04 


+0.05 


11.75 


LL 


630.249 


+0.02 


+0.04 


+0.06 


9.55 


LL 


633.084 


+0.01 


+0.02 


+0.04 


1.37 


LL 





0.019 
±0.011 


0.040 
±0.023 


0.061 
±0.034 






0* 


0.032 
±0.008 


0.070 
±0.014 


0.110 
±0.020 







Notes. Columns in the top part list the rest wavelength of the absorption lines and their abundance corrections for the three magnetic flux cases, 
the equivalent width obtained from the spectral synthesis based on our HD series, and the identifiers associated with each spectral line, following 
the convention introduced in Table [T] In the bottom part, () indicates the unweighted mean of the abundance corrections for the visible lines and 
corresponding rms for the different magnetic flux cases, and ()* refers to the values obtained when averaging for just the 5 lines in common with 
lAsnliindct al1l2000H) 
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line-to-line agreement found by those authors is not spoilt when 
the MHD approach is considered. 

Spectral lines in the IR: the average value of the abun- 
dance corrections for the IR iron lines we studied is indicated in 
Table[3](bottom row, marked with ()). We derive a negative aver- 
age correction in the range -0.02 < (Alog e(Fe) ) < -0.04 dex. 
However, the scatter around the mean is quite large for these 
lines, because of the scatter in their Lande factor (see Table [TJ 
and the relative importance that the Zeeman effect has for each 
of these IR lines. 



Table 3: Abundance corrections for the five IR spectral lines of Fei 
included in this study. 



A 
[nm] 


Alog 6(Fe) 

|dex] 
(50 G - HD) 


Alog e(Fe) Q 

|dex] 
(100G-HD) 


Alog e(Fe) Q 

[dexj 
(200 G - HD) 


W HD 
[pm] 


1558.826 


0.00 


0.00 


+0.01 


11.63 


1564.851 


-0.05 


-0.07 


-0.08 


9.57 


1565.287 


-0.03 


-0.05 


-0.06 


7.18 


1566.202 


-0.02 


-0.02 


-0.03 


19.88 


1566.524 


-0.01 


-0.01 


-0.02 


7.67 





-0.022 
±0.017 


-0.030 
±0.026 


-0.036 
±0.031 





5. Discussion 

First results using a 3D MHD approach to abundance analy - 
sis were presented in our previous paper dFabbian et al.lr20loh . 
where the study was limited to three iron lines. Here, we have 
aimed at extending those exploratory results using a much larger 
set of spectral lines, which allowed us to further improve the con- 
straints on the relevant derived effects on the solar abundance. 
However, it is important to note that, given the differential na- 
ture of our study, it is not our main aim to match the exact value 
of measured equivalent widths. Despite this, with the input solar 
iron abundance we adopted for the spectral synthesis and us- 
ing a sensible choice of magnetic flux, we find generally good 
agreement between our predicted MHD equivalent widths and 
those from observational data available in the literature for lines 
mark ed "A" in our Table [2] (e.g., the value by 
19951 for the 624.065 nm line and the values by 



Holweaer et al 



Blackwell et al 



1995 for the four remaining lines). In all cases we find that it 
is indeed appropriate to consider a magnetic flux in the range 
(B vert ) ~ 50 - 200 G at the corresponding depths of line forma- 
tion in the real Sun. 

Photospheric Fe lines in the solar case have the property that 
their strength decr eases with i ncreasing temperature in their for- 
mation layers (see[Grav 1992). Hence, the hotter average strati- 
fication caused by the presence of magnetic concentrations tends 
to decrease their predicted equivalent widths with respect to the 
HD case. The direct (i.e., Zeeman) effect acts in the opposite di- 
rection. In the visible range, the former effect dominates and sys- 
tematically smaller equivalent widths result. In the IR, Zeeman 
broadening is important and stronger lines (i.e., larger equivalent 
widths than in the HD case may result as a net effect. 

A dual picture thus clearly emerges from our results. For 
spectral lines in the visible wavelength range (Figs. [3)111 and[5]), 
Zeeman broadening is relatively small and the dominant influ- 
ence of magnetic fields on line formation is the line-weakening 
effect of the warmer average temperature structure in the revelant 
line-formation layers of MHD model atmospheres (see FigHJ. 
The 3D HD abundance results will need to be amended via pos- 
itive abundance corrections to increase their value to the one ap- 
propriate to the magnetic case considered most representative of 
solar conditions. On the other hand, Zeeman splitting becomes 
important in the IR region of the spectrum and tends to be able to 
more than compensate for the indirect line-weakening effect, so 
much so that for lines with particularly high Zeeman sensitivity 
(gL > 1.5, top panels of Fig. |5J, a net line-strengthening effect 
is noticeable. The 3D HD abundance will in this case need to be 
revised downward (applying negative corrections to it). 

T hese results are in li ne with our previous exploratory find- 
ings (Fabb ian et al.l 20101) A cross-check can be made since the 
three lines in that paper are included in the current study. The 
agreement is very good. Any small difference in the theoretical 



Notes. Same order as per columns 1-5 of Table [2] In the bottom row, 
we provide the (unweighted) mean of the abundance corrections and 
corresponding rms. 
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Fig. 7: Equivalent width percentage change of the 200 G case compared 
to HD, plotted against line rest wavelength. Filled diamonds are for 
the 28 iron lines listed in Table [TJ empty symbols for 6 lines of other 
elements which were used to complement our line list and for testing. 



equivalent width or in the derived abundance corrections (gener- 
ally in the sense of slightly more significant effects having been 
found now) are due to better spectral sampling and continuum 
coverage in the present study. We chose to decrease the number 
of snapshots and the spatial subsampling employed compared to 
our previous study, only after having made sure that this did not 
have a significant impact on the equivalent width of the lines. 

To confirm that the physical processes involved can explain 
the effects we found, we checked for any abundance correc- 
tion trends with the parameters characterizing the line formation, 
namely, rest wavelength of the line, lower level excitation poten- 
tial and oscillator strength of the given transition, level multiplic- 
ity (2s+l) and total angular momentum J, line centre formation 
height, Lande factor, and HD line strength. Here we present a fig- 
ure (Fig. [7} corresponding to the only trend retrievable, namely 
the sign of the abundance correction changes as one moves from 
wavelengths in the visible A range to the IR. In Fig.|7]we have ex- 
ceptionally included results for a few lines of elements other than 
iron. These are C 1 477.590 nm, 1 615.818 nm, Ni 1 676.777 nm 
(used in SOHO MDI magnetogram images), Mm 874.091 nm, 
Si i 1082.71 nm (this IR line is the strongest in the sample we 
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study, i.e., having an equivalent width > 500 mA), and finally, 
Mni 1526.27 nm (lAsensio Ramosll2009l :lA sensio Ra mos et alJ 
2009). We consider these additional lines of other elements be- 
cause, as shown in Fig. [71 they cover the region of intermediate 
wavelengths in the visible to near-IR where our sample has no 
iron lines, and provide hints that magnetic flux acts similarly also 
on lines of different elements than iron. Interestingly, the lines 
of the other elements in fact follow a similar trend to the one 
seen for iron lines. This would seem to suggest that the effects 
encountered can be generalized to other elements according to 
considerations on a given line's formation height and tempera- 
ture sensitivity, on the one hand, and Lande factor and wave- 
length, on the other. We also see that the percentage change in 
equivalent width decreases smoothly from the lower to the upper 
wavelengths range limits for those lines. 

In general, we find significant differences in the effects for 
different lines. Correspondingly, the abundance corrections vary 
from quite small (|Alog e(Fe) | g 0.02 dex) to significant (up to 
|Aloge(Fe) | ~ 0.15 dex in the 200 G case) for all but a few 
of the spectral lines analysed. Lines affected to a level below 
0.02 dex represent only rare exceptions. 

We focussed also on the subset of lines that we have in com- 
mon with the absolute solar iron abundance determination of 
lAsplund et al.l (2000b|) These all have gi > I. The abundance 
corrections found here for them are among the most significant 
of those we find to be applicable to the lines we have studied. 
Interestingly, for this subset of lines, the effects due to the pres- 
ence of magnetic fields can be seen as linked to their formation 
height. That is, they follow the trend seen in Fig. [T| according 
to which lines formed between a height of 200 and 400 km feel 
increasingly strong indirect effects. The abundance corrections 
for lines in this subset follow the trend expected, based on their 
formation heights (see values in iG urtovenko & Kostikl I1989L 
Table 2). Oth er spectral features used in the study bv Asplund 
et al. (2000b) are also likely to be affected by significant abun- 
dance corrections, in particular for lines with a small Lande fac- 
tor, for which the indirect effect found here can act undisturbed. 
Thus, the last word on the solar iron photospheric abundance has 
probably not yet been spoken. 

Our results also urge more investigation into the magnitude 
of the magnetic effects for other chemical elements than iron. In 
particular, this is true for those elements whose abundances have 
in the last decade undergone a revision that has caused large dis- 
agreement between standard solar models and helioseismology 
constraints. It is crucial that also oxygen, carbon, neon, nitrogen 
and silicon - i.e. those elements playing a major role in deter- 
mining the total solar metallicity Z/X - be redetermined via tar- 
geted studies simultaneously including the most accurate atomic 
and observational data and non-LTE, 3D, and magnetic effects. 

Since all spectral lines have a more or less pronounced sen- 
sitivity to the temperature values of the layers where they are 
formed, and since direct magnetic sensitivity of the lines, if 
present, induces line broadening, one can make a broad division 
into the following expected general behaviour, depending on the 
chemical element and ionization stage considered: 

- Spectral lines for which magnetic broadening acts in the op- 
posite sense to the indirect temperature effect. This is the 
case of both neutral and singly-ionized iron lines in the Sun 
and of other elements having the same temperature sensi- 
tivity. In this category, there will be some rare instances 
of close-to-negligible net abundance corrections. In general, 
though, the direct and the indirect effect will not balance out 
exactly. The net resulting abundance correction can be ei- 



ther positive or negative, with a sign inversion expected when 
moving from UV to IR wavelengths, i.e., towards the spec- 
tral range where Zeeman broadening becomes important. 
- Spectral lines for which the indirect temperature effect acts 
in the same direction as magnetic broadening. These will 
tend to become stronger in MHD simulations, thus the rel- 
evant abundance corrections to be applied to 3D HD abun- 
dance determinations will always be negative. Still, their 
magnitude should be investigated on a case-by-case basis. 
The net abundance corrections will be maximized, because 
of the direct and indirect effects working together to rein- 
force the spectral line. Although Zeeman splitting for most 
lines of chemical elements other than iron may not be as im- 
portant, the spectral features in this category will still tend 
to become stronger in the warmer MHD temperature strati- 
fication thanks to the indirect effect. One important example 
are lines of neutral oxygen (since that chemical element re- 
mains mostly neutral in the solar photosphere due to its much 
higher ionization potential compared to e. g. iron), for which 
our preliminary abundance results of work included in an 
upcoming publication seem to indicate significant effects. 

Given the differential nature of our investigation, we expect 
that the neglect of non-LTE effects on Fei lines should not 
change our conclusions significantly. In fact, if anything, non- 
LTE corrections for neutral iron are expected to push the so- 
lar abundance determination even higher due to over-ionization 
caused by the near-UV radiation field, hence producing under- 
population of the neutral iron atomic levels and thus weaker 
emerg ent line profiles compared to LTE due to line opaci ty 
deficit (IColletet al.l ( 20051) : IShchukina & Truiillo Buenol ( 200 If) ). 
In fact, attempts at improving the 3D HD LTE-based Asplund 
et al. (2000b) study also in this sense have already pushed the 
solar iron abundance back up to log(Fe) Q = 7.50, while achiev- 
ing better agreement with helioseismology, and while preserv- 
ing the agreement between values based on lines from neutral 
and singly-ionized Fe and the small line-to-line scatter and ab- 
se nce of trends with excitation potential (see Table 1 and Fig. 6 
in lAsplund et al. 20091). Based on that result and on our current 
findings regarding the systematic abundance corrections applica- 
ble due to magnetic fields (with likely additional e ffects expected 
in 3D non-LTE, see iHolzreuter & Solan ki 2012), the real pho- 
tospheric solar iron content may well be log(Fe) e > 7.50, i.e. in 
a relatively high range compared to the best estimates obtained 
in the past decade. A less "anemic" Sun would better agree with 
new, accurate deter minations of the solar Si cont ent accounting 
for non-LTE effects ( Shchukina & Sukhorukovl 2012f) 

Both of these changes, together with hints that 3D-based re- 
visions in other chemical elements may seemingly not need to be 
drastic compared t o the ID case (e.g.. see th e results qf^Caffau 
et al. 2009 based on CO s BOLD 3D granulation models), should 
help alleviate the discomfort that helioseismology has been feel- 
ing towards low values of the solar metal content. In any case, 
it will be very interesting to investigate how 3D, non-LTE, and 
magnetic effects interact, and whether there may be some ampli- 
fication or reduction in their magnitude when looking at a fully 
consistent picture in order to carry out an absolute abundance 
determination. 

A further remaining uncertainty revolves around choosing 
the MHD model that is thought to best represent the real Sun, in 
order to estimate the most appropriate abundance correction to 
be applied to the abundance derived in the HD case. In particular, 
further testing of the realism of the MHD models must be carried 
out in terms of geometric configuration and depth profile of the 
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magnetic field. Our non-magnetic model has significantly lower 
values of the photospheric temperature for the average stratifica- 
tion at log T5oonm < compared with the temperature structure 
of the s emi-empirical solar atmosphere derived bv Holweger & 
Mueller ( 1974) For those same optical depths, our MHD mod- 
els are progressively closer to it as magnetic flux is increased, 
while very good agreement is generally re ached for deep lay- 
ers. Compared to the theoretical results of IStein & Nordlundl 
(1998) (s ee their Fig. 15. also shown in Fig. 14 of Nordlund 
et al. 20091), our HD atmospheric temperatures achieve similar 
to significantly better agreement (in photo spheric and deep lay- 
ers, respectively) with the stratification of iHolweger & Mueller! 
( 1974) A number of observational results are certainly helping 
to constrain the possibilities, and comparisons with future, high- 
quality, targeted solar data are set to improve this scenario even 
further. In this sense, we have embarked on a detailed compar- 
ison and testing in terms of line parameters and thermodynam- 
ics, of centre-to-limb variation (CLV) and Stokes parameters, as 
well as of line fitting and absolute abundance derivation, but this 
goes beyond the scope of the present paper and is part of work 
in progress. As an example, the tests we carried out and dis- 
cuss in a separate article (Beck et al. 2012, submitted), show 
that the agreement with the observational data we studied of 
predicted properties (for spatially- and spectrally-degraded spec- 
tral lines), such as the rms fluctuations at different line depres- 
sion levels of bisector position/velocity, bisector intensity, and 
line width, may at least in some cases improve when including 
magnetic fields in the simulations. It will also be important to 
continue along the line traced by initial effort s of using Stokes 
spectra from magneto-convection simulations dStein et al.l 201 1 ) 
to make comparison possible of theoretical predictions with the 
best available current and future spectropolarimetric observa- 
tional data. 



6. Conclusions 

We have demonstrated that appreciable effects on a set of iron 
spectral lines are found when introducing magnetic fields in 
3D, time-dependent radiation-hydrodynamics simulations of so- 
lar surface convection. Even the most advanced chemical com- 
position studies so far, using 3D convection modelling, may still 
be affected in a significant way in terms of metallicity determi- 
nations owing to their non-magnetic assumption. We here deter- 
mined the associated abundance corrections to be applied when 
abundances are derived via a purely HD approach. In the follow- 
ing, we briefly highlight our main findings. 

We have presented evidence that the magnetic field causes 
both Zeeman broadening in the lines with non-zero Lande fac- 
tors and also changes in the average temperature stratification of 
the cells. Since the majority of the lines have non-zero Lande 
factors and are also temperature-sensitive to some extent, no 
spectral line is in principle exempt from relevant abundance cor- 
rections. The effects are therefore ubiquitous. The direct (line- 
strengthening) effect of magnetic field on line formation via 
Zeeman broadening is generally (except in the IR) less important 
in magnitude than the indirect effect of temperature stratification 
changes. 

The abundance corrections are approximately proportional 
to the average unsigned magnetic field strength (|B velt |). For 
(Evert) = 100 G, the abundance corrections are |Alog(,Fe)| ~ 
0.04 dex for our full selection of iron spectral lines in the visible 
wavelength ra nge and |Alog(Fe)| ~ 0.07 dex for the subgroup in 
common with lAsplund et al.l (2000b) 



If, as a number of recent results seem to suggest, the "quiet" 
solar photosphere is in fact teeming with hidden magnetic flux, 
the abundance of iron based on HD modelling will need to be 
revised upwards. It is likely that a more accurate estimate for the 
solar iron content may be log(Fe) > 7.50. 
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